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ABSTRACT - We prove that there are no Petrov type III space-times on which the 
conformally invariant (self-adjoint) scalar wave equation or the non-sclf-adjoint scalar 
wave equation satisfies Huygens' principle. 



RESUME - Nous prouvons qu'il n'existe aucun espace-temps de type III de Petrov 
sur lequel l'equation invariante conforme des ondes scalaires ou l'equation des ondes 
scalaires non-auto-adjoint satisfait au principe de Huygens. 



1 Introduction 



This paper is devoted to the the solution of Hadamard's problem on Petrov type III 
space-times, for the conformally invariant scalar wave equation 

Du + lRu = 0, (1) 
o 

and the non- self-adjoint scalar wave equation 

Du + A a d a u + Cu = 0. (2) 

^Published on Ann. Inst. Henri Poincare (A) Phys. Theorique, 71, 595 (1999). 
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In the above equations □ denotes the Laplace-Beltrami operator corresponding to 
the metric g a b of the background space-time V 4 , u the unknown function, R the Ricci 
scalar, A a the components of a given contravariant vector field and C a given scalar 
function. The background manifold, metric tensor, vector field and scalar function 
are assumed to be C°° . All considerations of this paper are entirely local. 

The homogeneous equations (0 and (J2J) satisfy Huygens ' principle in the sense of 
Hadamard [H] if u(x) depends only on the Cauchy data in an arbitrarily small neigh- 
borhood of the intersection between the backward characteristic conoid C~(x) with 
the vertex at x and the initial surface S, for arbitrary Cauchy data on S, arbitrary 
S, and for all points x in the future of S. Hadamard's problem for (0) and (J2J) is that 
of determining all space-times for which Huygens' principle is valid. We recall that 
two equations (0) are said to be equivalent if and only if one may be transformed into 
the other by any combination of the following trivial transformations: 

a) a general coordinate transformation; 

b) multiplication of the equation by the function exp{— 20(x)), which induces a con- 
formal transformation of the metric 

9ab = e 2(p g ab ; 

c) substitution of Xu for the unknown function u, where A is a non- vanishing function 
on V4. 

We note that the Huygens' character of (J2J) is preserved by any trivial transformation. 
In the case of (0 the trivial transformations reduce to conformal transformations with 
A = e*. 

Carminati and McLenaghan [21] have outlined a program for the solution of 
Hadamard's problem for the scalar wave equation, Weyls' neutrino equation and 
Maxwell's equations based on the conformally invariant Petrov classification of the 
Weyl conformal curvature tensor. This involves the consideration of five disjoint cases 
which exhausts all the possibilities for non-conformally flat space-times. Hadamard's 
problem for (JTJ) and © has been completely solved for Petrov type N space-times by 
Carminati and McLenaghan |5| and McLenaghan and Walton [T^]. Their results 
may be summarized as follows: 

Any non- self- adjoint equation on any Petrov type N background space- 
time satisfies Huygens ' principle if and only if it is equivalent to the wave 
equation □« = on an exact plane wave space-time with metric 

ds2 = 2dv{du + [D(v)z2 + D(v)z 2 + e(v)zz]dv} - 2dzdz . (3) 

For Petrov type D space-times the following result was obtained by Carminati and 
McLenaghan McLenaghan and Williams [201 and Wiinsch |26j : 
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There exist no Petrov type D space-times on which the conformally invari- 
ant scalar wave equation (QJ) satisfies Huy gens' principle. 



In the present paper we complete this program for the conformally invariant scalar 
wave equation © on Petrov type III space-times by proving the following theorem: 

Theorem 1 There exist no Petrov type III space-times on which the conformally 
invariant scalar wave equation (QJ) satisfies Huygens' principle. 

The results on type N and type D space-times described above and Theorem 1 lend 
weight to the conjecture which states that every space-time on which the conformally 
invariant scalar wave equation satisfies Huygens' principle is conformally related to 
the plane wave space-time © or is conformally flat |H |3] . 

Hadamard's problem for the general non-self-adjoint equation (j2J) may now be 
solved with the help of Theorem ^ and the results of Anderson, McLenaghan and 
Sasse PP where the following theorem is proved: 

Theorem 2 Any non-self-adjoint scalar wave equation |^) which satisfies Huygens' 
principle on any Petrov type III background space-time is equivalent to the conformally 
invariant scalar wave equation (QJ). 

Combining these two theorems we obtain 

Theorem 3 There exist no Petrov type III space-times on which the non-self- adjoint 
scalar wave equation satisfies Huygens' principle. 

The corresponding problem for the Weyl neutrino equation and Maxwell's equations 
is solved in [T%] . 

The starting point of our proof of Theorem Q is the paper by Carminati and 
McLenaghan ,7 a , where the following results are obtained for Petrov type III space- 
times: 

Theorem 4 The validity of Huygens ' principle for the conformally invariant scalar 
wave equation on any Petrov type III space-time implies that the space-time is 
conformally related to one in which every repeated principal spinor field oa of the 
Weyl spinor is recurrent, that is 



'A;BB 



BB i 



(4) 



where I BB is a 2-spinor, and 



ABCD;EE 




0. 



(5) 




0. 



(6) 
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Theorem 5 If any one of the following three conditions 

*ABCD;EE^ B L D L E O E = 0, (7) 

*ABCD i Bi:L A i B o D <W t = 0, (8) 

Vabcd;EeL A l B l D o E o e = 0, (9) 

is satisfied, then there exist no Petrov type III space-times on which the conformally 
invariant scalar wave equation |7J) satisfies Huygens' principle. 

It is important to note that these earlier results solve Hadamard's problem under 
what have proved to be fairly strong assumptions (namely, that one of (J2J), (jEJ), or 
dlU is satisfied). The purpose of the present paper is to make the analysis completely 
general by removing these assumptions. We follow the conventions of [7j, and use the 
results established there to obtain (most of) the basic equations needed for the proof 
of Theorem 

In Section El we give the necessary conditions for the validity of Huygens' principle 
that will be used in this paper, and give a brief summary of their implications. From 
these necessary conditions, we derive the further side relations needed for our anal- 
ysis in terms of the Newman-Penrose scalars. The key to our proof is the six-index 
necessary condition obtained by Rinke and Wiinsch [22] which was not used in [7j. 
In Section |21 we examine these side relations in the case $xi — and show that they 
lead to a contradiction. The proof of Theorem Q is completed in Section where the 
case $n 7^ is treated. 

It is worth mentioning that the tools of computer algebra are used throughout 
this paper. Initially we employ Maple package NPspinor [5J GUI to extract dyad 
components of spinor versions of the necessary conditions, and then to manipulate 
the resulting expressions in Newman-Penrose form. In the case $n = we use 
the Grobner basis package of the Maple system to explicitly determine solutions of 
systems of algebraic equations. Finally, for the case $n ^ we use the GB package 
of Faugere ^T] to examine the solvability of a somewhat larger system of algebraic 
equations. 

2 Formalism and Basic Equations 

The necessary conditions for the validity of Huygens' principle for (0) which we employ 
are given by 

(III) S abk . k - ^C k ab l L kl = 0, (10) 

(V) TS (3C k ab l . m C kcdl + 8C ab - c Skid + 40S' a fe S c dk — 8C ab S k ic-,d 
— 24C ab S c dk;i + 4C ab Ci n + l2C k ab l C™ 

cdl^-'km 

= > (11) 
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(VII) ts(q 

abcdef abcdef + abcdef + abcdef + Q 



(5) \ 



(12) 



where 



Q 



(i) 

abcdef 



3C k ' n h ' m c Ckdel;mf + C^ab ;cd(105fc/ e; / + 6S e fk;l) + ^^abk;cSde 
cdk;efLlm + 5Ckcdl;meL n f + 7C" 

cdkde-Lmf 



ab ; c*- 



-C\ b l (3C 

+135 , M C;( iI/ e j- + 12S c dk-,iL e f + 71S 1 , 



cdk;eLlf) , 



(13) 



^abcdef 



C k a b ;c {Skld;ef + 3Sdek;lf + 2Sabk;cdSef k ~ 5S a bkS c dL e f 
1 



^'-C'ab ;c (2C m kld;eL m f + 3(7™ dek;lL m f + SkldLef 



+3Ckde m ;fLl m + ^SdekSlf) ~ C* a b \Cj 
+ S c dkL(le-J) — —R ]C Ckdel;f) , 



kcd ;e-^ J (lm;f) 



(14) 



^abcdef 



/-ik I (rifi mn /~i -t n/imn >h( i 9(1^1 mo N 

^ 06 c^lnmd;ef ~ lu <^ c dykefl;mn "T ^UL// C rf Okme;f ) 

cClmndI>ef 



"O^fc a^lmnb^ cd ;ef T O a fe l,' ( -'fc tA^lrrmd-L', 



— lOCkeflC 



m n r 



(15) 



^abcdef 



— <^ afe c;d^lmnd;ef + U^lcd ;e*->fcm/ + '^Icd ;k>->efm 

76 404 

— — t-'cfcZ ;dJefm — ^cdk^efl ) + OL'fc a<-^ bc;d^ lefra;n ■, UOJ 



^abcdef 



/~ik m t t 

— ^ ab Wed ■L/hm.-L'. 



km^ef 



where 



(17) 



'S'afec := -^o[6;c] > 

-^ab := — -Rab + -^9ab ■ 




(18) 
(19) 

(20) 



Here A a := g a bA b , R a bcd denotes the Riemann tensor, C a b c d the Weyl tensor, R ab : = 
g cd R ca bd, the Ricci tensor, and R := g ab R a b the Ricci scalar associated to the metric 
g a b- The conditions III, V and VII are necessarily conformally invariant. Spinor 
versions of conditions III and V, and the conventions used for conversion from the 
original tensor form, are given in [7j. 

Mathisson JE], Hadamard and Asgeirson [2] obtained condition III for (J2J) 
in the case g iJ constant. Condition /// was obtained in the general case for (0) by 
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Giinther [T3]. Condition V was obtained by McLenaghan [T7j in the case R a b = 0, 
and by Wiinsch [23] for the general case. Condition VII was obtained by Rinke and 
Wiinsch [22]. 

Petrov type III space-times are characterized by the existence of a spinor field o A 
satisfying 

^abcdo c o d = , V A bcdo d ± 0. (21) 

Such a spinor field is called a repeated principal spinor of the Weyl spinor and is 
determined by the latter up to an arbitrary variable complex factor. Let i A be any 
spinor field satisfying 

o a l a = 1 . (22) 

The ordered set oa, i-a, called a dyad, defines a basis for the 1-spinor fields on V4. 

It was shown in [Zj that the necessary conditions III and V imply that there 
exists a dyad {0,4, l a } and a conformal transformation <ft such that 

K, = a = p = T = e = 0, (23) 

^0 = *i = *2 = ^4 = , * 3 = -1 , (24) 

$00 = $01 = $02 = A = 0, (25) 

Da = D(3 = Dtt = 0, (26) 

«J$ n = D$ n = 0. (27) 

We notice that the expressions (|24|) determine the tetrad uniquely. On the other 
hand, conditions (|2~H|) are invariant under any conformal transformation satisfying 

D0 = 0, 50 = 0, (28) 

which implies that we still have some conformal freedom. Under a conformal trans- 
formation we have |24j : 

$ n = e-**$ n - (29) 

Thus, we can choose <fi such that 

$11 = c, (30) 

where c is a constant. The conditions (j28|) are satisfied in view of ([27]) . 

Let us now derive some side relations that follow from the previously obtained 
equations (|2*3|) - (p?0|) and the necessary conditions (fTU|) - ([T^jl ; these will be required 
in the analysis of the following sections. We may assume that a/3ir 7^ 0, since the 
case in which this is not true was already considered in [7j. By contracting condition 
III with l a o b T ab we get 

5/3 = -/3(ju + 0). (31) 
Using the Bianchi identities and the above conditions, we obtain 

D$ 12 = 27f$ n , (32) 
D$ 22 = -2(P + p) + 2$ 21 7T + 2$ 12 7T , (33) 
5<t> 12 = 2a + AW + 2A$n - 2a$ 12 , (34) 
5<S> 12 = -2(3 + 2p$u - 2^$i 2 . (35) 
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The Ricci identities provide the following relevant Pfaffians: 



Z>7 = 


--Wa + (3n + $11 , 




(36) 


DA = 


: (l/2)8a-(ll/2)0a- 


f 7T 2 - 2na - lln/3 - (3/2)a 2 , 


(37) 


5w = 


DX — W 2 — Wa + 7T/3 , 




(38) 


DT7 = 


: AW + Wp, + X7C + IFj 


- 7T7 - 1 + $12 , 


(39) 


5a = 


6/3 + aa + 00 - 2/3a 


+ $11, 


(40) 


5tt = 


Dfi — Wn + na — (37T . 




(41) 



We can obtain useful integrability conditions for the above Pfaffians, by using Newman- 
Penrose (NP) commutation relations. By substituting them in the commutator ex- 
pression [8, -D]$22 — [A,-D]$i2, we get 

50 = -2$n -0a- 40rf - 2DJ1 -0/3 + 2txtx. (42) 

By contracting condition V with l abcd I ab o od , we find 

20071 + 120a + 6rra + 3a 2 + 5a + 25tt + 10 + 2 = . (43) 

By substituting (f3T)|) into this equation we get 

6(2W + a) = -20n(3 - 110a - 6Wa - 3a 2 . (44) 

From (gOI), flUJ) and (jSH we then obtain: 

5(2% + a) = 2na + aa - Q/3n - 3(3a - $n . (45) 

By contracting condition V with i ABC o D I ABC o D , we find 

—6 5 7T — 15aW — 10 aa — 68tcW — 15 7r 57 — 35 a — 126 (3 (3 

+5D7 + 10 Dp - 240a- 35a -65W -15 50 - 30W 

+5D 7 +IOD/1 - 155/3 -24/3a-3/37r-4$n = 0. (46) 

Using (j37)|). (j4*T|) and the complex conjugate of (jI2J), we get 

9$ u + lO/ftf + 5[D/i + D/7] - 2vra - 1200 - 2aa - 16Wn + 10(3tt - 2na = . (47) 

On the other hand, the NP commutator [6, 5}(a + 2n) = (a — (3)5(a + 2tt) + (—a + 
f3)5(a + 2ir), yields the following expression 

2ir(30 + 22n0a + A3tt0W - 22Wtt 2 + 0T>n + 22ttDji + 12a0a 

+60<1> U - 12aWn + lla$n + 187r$ n + 24w0a + 12aD/J = . (48) 
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Eliminating D/i between (j4T)) and P8|). and solving for D/i, we get 

Du = — f 108 vr$ n - 44 vr 2 a - 24 vr a 2 - 68 vraa 
5 v 

144 app + 53 a$n - 274 nPp + 120 n(3a - 24 a 2 a 

-242vfvr 2 + 220/?vr 2 - 176 aW p - 60 a/3a- 30/?$u 

+5vr^- llOavr/?) /(-/3+ 12a + 22vr), (49) 

where we have assumed that the denominator of the expression above, given by 

d x := -J3+ 12 a + 22n, (50) 

is non-zero. The case d\ = will be considered later. 

Substituting expression (j4*9~|) for D/i into (|4T|) we obtain 

D/Z = — - (90tt$ii + 12/3/3 2 - 10/fvf + 55a$ n + 122apW - llOvfvr 2 - 60avrvr 
5 v 

+62apa + 21/3$ n + 231vr/?7r + 112371$) /(-p + 12 a + 22 vr) . (51) 

One side relation can now be obtained by subtracting the complex conjugate of (f4"9|) 
from 1)51 j) . We obtain: 

1 — —9 

Si := -(720a 2 / 3a + 2904vr 2 vr 2 - 12/3 2 /3 + 288 a 2 a 2 + 528 vr 2 a 2 + 528 vr 2 a 2 
5 

+2420 a vr p vf + 3056 a vr /3 /3 + 2888 a a vf vr + 1320 a vr /3 a + 1320 vf (3 a 2 
+ 1606 a ppa + 5802 Wit pp + 1320 fa^a + 2420 W/3 a vr + 1320 vr/?3 2 
+3056 vr /5 a p + 2552 vr vr 2 a + 305 /? $ n a + 570 /3 $ n vr - 51 /3 $ n /5 
+24 a $11 a - 86a$n vr + 305 3$n/? + 816 vraa 2 + 2552 a vf vr 2 
+816 avra 2 - 86vra$ n - 396vfvr$ n + 720 paa 2 + 570 7f/3$n/ 

((-12a-22vr + /3)(12a + 22vf- / 3)) = 0. (52) 

We notice that ()49)1 and (|HT|) have the same denominator. Thus, if we keep these 
expressions for D/i and D/Z, the Pfaffians SP, 5a, 5tt, given by (JI2J), (J3%|) and (JHJ, 
respectively, and their complex conjugates, also have the same denominator. This 
procedure is crucial to keep the expressions to be obtained from the integrability 
conditions within a reasonable size. Except for 5a, all Pfaffians involving 5, 5, 
applied to a , j3 , vr are explicitly determined. 

The following expression for 5a can be obtained from the NP commutator [5, 5]P = 
(J[- /i)D/3 + (a-P)5p+ (P -a)Jp : 

5a = (2vrD/Z - 2S(Dp) - 3vra2 - 8/3D/Z - Ua]k - 2vfvr 2 

-4avrvr - 4/3$u - UTpW) /vr . (53) 
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By substituting (}4*9"j) and (J5TJ) into this equation we get: 

5a := - -(19519 $n ~f a + 8570 $n /? a 2 + 1950 $ u vr a 2 + 35850 $n vr /3~ 2 
5 

+2900 $ n vr 2 /3 - 3180 vr 2 /3~ 2 vr - 210 a 3 ~j3 tt + 150 a 3 ~j3 a - 1307 a 2 ~f W 
-180a 2 p p 2 + 628a 2 /3~ 2 a - 1280 7r 2 /3 2 a + 1950 a/3~ 3 a + 4668 a/3 3 /3 
+3520 /f a 7T + 8160 P tt /f + 330 p 3 W tt + 975 $ u a 3 - 860 /3~ 3 $n 
+ 17950 $n /? a vr - 420 a 2 p W tt + 300 a 2 P tt a - 41 16 a /f W tt 
+588 an a p 2 + 175 p 3 an) /(p 2 W - UaPw + 10 P a a + 20 p tt a 

+65 a $n - 12 /3 p 2 + 130 tt $n - 31 P $ u - 28 /3 tt tt - 6 /? 2 a) , 
where, for now, we assume that the denominator in the expression above: 

d 2 := p 2 W- UapW + 10 /3a a + 20/3 vr a + 65 a $n - 12/3 /f 
+ 130tt$ii -31/? $n - 28/37f7r-6/fa, 

is non-zero. 

Contracting condition VII with l ab c de f o bcd t aef yields 

VTiia := 648/3 2 avr - 165 vr tt 5/3 - 40 tt D$ 2i + 1557 5/3 a tt 
+270 a /3 5vr + 144 5a a /3 + 378 a /3 D7 + 2058 a 2 vf vr 
+570 D/i 5/3 + 630 D/i /f + 21 DA 5a - 201 D7 5vr 
-36 vr 2 5a - 954 5/3 p 2 - 864 5/3 5^ - 864/3 p 3 - 132 5a 5a 
+ 1116a 3 a- 366a 2 5a- 900a/35/3 - 20aD$ 2 i + 1454 $ U 7r/3 
+771 5vf vr^ + 423 5a a tt + 939 a a 5tt - 1038D7Ja7r 
-90AD(D7i) +45 vr 5(5/3) - 15 tt 5( D7 ) - 30 vr 5( D71) 
+84 a p5w - 774 5a /3vr - 570 5/3 a tt + 180 5a /3 /3 
+ 135 D( 5/3) A — 45 A D( D7) - 186 5/3 a 7f + 570 a /? D/I 
+30 DA a 7f — 39 DA a 7r - 78DA7f7r - 744aa5/3 - 924aa 2 /? 
-387 tt 2 a P - 21 DA a a + 351 D7 tt /3 - 1773 /3 vr 5tt 
+216avr57r + 870 D/ia/3 - 198 5a a/3 + 762 5aa a 
+378 p a 5/3 + 828 a 2 /3/3 - 1674 a p 5/3 - 594 /3 /3 5/3 + 636 vr 2 a a 
+2289 a 2 tt a + 36 DA a /3 - 660 a vr 5a + 1773 /f 7f tt - 555 D7 a tt 
+2187 /3 tt /3 2 + 2619 p tt Jp - 888 a p 2 a - 126 5a a p 
+3681 p tt a P - 2055 a tt a /3 - 6951 /3 7f a tt - 1476 a /f /3 
+ 150 5a tttt - 63 DA /3vr - 30 /f a vr + 934 $ n a tt - 465 tt 2 /3tt 
-1608 a 7r57f + 63 DA/3/3 + 20 /3D$ 2 i + 567vr 2 pp 
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-426 5n5n + 18 DA D/Z - 39 vr 2 D7 + 39 vr 3 a - 78 vr 2 D/Z 

+78 tt 3 vf + 9 DA D7 + 324 5a 55 + 324 5a ~f 

-189 5a 5tt + 276 50 D 7 + 306 ^ 2 D 7 + 432 a 3 ¥ 

+1080 a 2 50 - 144 5a 5vf + 360 5a 50 - 372 a 2 5W 

+630 5/3 5tc + 630 /3~ 2 5tt + 528 D/Z tt + 408 tt tt 5tt 

+246 a W 5tt + 639 5a tt + 156 a D7 - 2592 a 2 /3 ¥ 

+204 5a a W + 117 /3 /3 5n - 657 5/3 tt /3 + 1 128 tt 2 a W 

- 108 5a D7 - 240 5a D71 - 720 a 2 D/Z - 324 a 2 D7 

+567 50 5tt - 30 DA 5tt - 27 DA 50 + 1 17 vr 2 5^ 

-378 DZZ57T - 102 $ n DA + 182 $ n tt 2 + 366 $ n 5tt 

-600 $11 5/? - 600 $11 /f + 1008 $ n a 2 + 336 $ n 5a + 1704 $ n a /3 

+30ttD(57) - 90 tt 5(55) + 60 7rD(5zz) -90D(5^)// 

-90ttD(A^) - 660 ^ vr D 7 - 270 a 5(5^) -30a5(5a) 

-60a5(5vf) +30^5( D 7 ) + QO05( Dfi) + 90 a5( D7) 

+180a5( D/Z) - 180 /3 5( 5?r ) - 90/3 5( 5a) + 45 tt5( 5a) 

+90tt5(5¥) + 30^D(5 7 ) + 60^D(5/i) + 15AD(5a) 

+30AD(5¥) + 60a5( 5a) -90 05(50) -90^D(Aa) 

-180^D( Att) + 180/35(5^) + 120 a 5( 5tt ) + 60 7 D( 5a ) 

+1207D(5tt) - 18Oa06>ir-6O0Dpir-315 050ir 

+180/357T7T = 0. (56) 

The second-order terms D(5 7 ), D(5 7 ), D(5/Z), D(A/3), D(Aa), D(A7r), can be ex- 
pressed in terms of known Pfamans and 5a, by using the NP commutation relations 
involving each pair of operators. After the substitutions we obtain 

VII 13 := % (760 2 vr n 5a + 9600 0<S> U 5a tt + 25040/3 5a a tt 2 
5 

-23440 0W7r 2 5a + 2400 vr 5a 2 - 240 2 5a a tt 

+2880 a 2 5a + 8600 a $ n 5a tt + 360 a /f W 5a 

+600 a 5a a /f + 24880 a 5a a n - 221Q0 a 0Wtt 5a 

+8940 a $ n 5a 0- 5280 a 2 /3vf 5a + 6240 ^aa 2 fo- 9885300 a $n tt 2 /3 

+ 113040 a00 A - 1200 $ n 5a a 2 + 48400 a tt 3 $ n + 40200 $ n tt a 3 

+22400 $ n 5a tt 2 - 720 ^ 3 /3 5a - 300 5a a ^ 3 + 142635 a0 3 $ n 

-910512 a 2 3 - 351912 2 a 3 tt - 5085800 ^ vr 3 $ n - 1467420 $ n a 3 

-1760 $ n 5a 2 - 3600 $ n a 4 + 35640 3 a 2 W- 360060 a0 3 a 2 

+ 120000 $ n tt 2 a 2 - 633396 $ n /f a 2 + 8640 a 3 /3/f - 15840 a 4 /3> 

-155352 a 3 /f a + 18720 a 4 /3 a + 176880 tt 2 3 W + 205920 0tt0 A 
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+85800 jf a vr + 47100 a ^ 4 a + 227170 $ n tt /f - 932620 $n vr 2 /f 

-2702400 tt 2 /3/5 3 - 1041760 vr 3 /f a - 2977440 vr 3 /f ¥ - 1064800 vr 2 /f a 

-835840 a 2 vr a /f - 6584360 $n /? a 2 n + 74640 a 3 p vr a 

-66480 a 3 P W vr + 75120 vr 2 /? a 2 a - 70320 vr 2 a 2 /3 vr - 2193544 a 2 /f vr vr 

+7200 a 2 (3 Ti P 2 - 4461264 a vr 2 /f vr - 3128928 a /? vr /f 

- 1237920 a /f a vr + 161720 a jfTir - 1662504 a $u vr/f 

-1572528a7r 2 /3 2 a-25/3 4 $ n )/(-/3 + 12a + 22 7r) 2 = 0. 

Solving this equation for 5a we get 

5a := - (227170 tt^^h +40200 vr a 3 $n + 176880 vr 2 /f W + 18720 /3a 4 a 
-1467420/3 a 3 $n - 5085800 /3 vr 3 $n + 85800 n]f a - 1064800 vr 2 /f a 
+8640/f a 3 /3- 932620 7r 2 /f$n - 1 5840/3 a 4 vf + 205920 vr/? 4 /? 
-2977440 7r 3 /3 2 7r- 1041760 vr 3 a /f - 2702400 vr 2 /3/f - 360060 a 2 /f a 
-351912 a 3 /f W - 155352 a 3 /f a - 633396 a 2 /f $n - 910512 a 2 
+35640 a 2 /f vr + 47100 /faa + 142635 /f a $n + 113040 /f a /3 
-3600 a 4 $n - 1572528 tt 2 /fa a - 1662504 Tip 2 a $ n + 161720 Tr/favr 
-835840 vr p 2 a 2 a - 1237920 np 3 aa- 4461264 tt 2 a /f ¥ 
-2193544 a 2 irp 2 W+ 7200 ]f a 2 vr /3 - 3128928 /3~ 3 a vr p 
-9885300 Ptc 2 a $n + 75120 /3V a 2 a - 70320/3V vra 2 
-6584360 P a 2 vr $n + 74640 /3 a 3 vr a - 66480 p a 3 vf vr + 48400 vr 3 a $n 
+ 120000 a 2 vr 2 $n -25^ 4 $n) /(-1760^ 2 $n - 300 ^ 3 a - 720 /3^ 3 

+360 /f a vf - 240 /f a vr + 2400 p vr /f + 600 a /f a + 2880 a /f /3 
+760 p 2 W vr - 5280 a 2 /3 ¥ - 22160 pw a tt + 25040 tt 2 a /3 
+8940 /3a$n + 9600/3vr$n + 6240 a a 2 /3 - 23440 vr 2 /3 vf 
+24880 a vr a /3 + 22400 vr 2 $n - 1200 a 2 $n + 8600 tt a $u) , 

where the denominator of ([58)1 . 

d 3 := -1760 ^ 2 $n - 300^ 3 a- 720/?^ 3 + 2880a^ 2 /3 
+360 /3~ 2 a vr - 240 /3~ 2 a vr + 2400 /3 vr /f + 600 a /f a 
+760/3 2 7f7r - 5280a 2 /3¥- 22160/3>a7r + 250407r 2 a/3 
+8940 p a $n + 9600 /3 vr $ n + 6240 a a 2 /3 - 23440 tt 2 /3 ¥ 
+24880 a vra^+ 22400 vr 2 $n - 1200 a 2 $n + 8600 vr a $n , 

is assumed to be non-zero for now. 
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Subtracting ()54j) from (|58|l and taking the numerator, 

Ni : = 684288 /f p 2 a + 286000 tt 3 a $^ 

-7290900 tt 2 /f $^ - 7913100 a 3 $^ /? + 165600 /f a 2 a 

-655680 tt 3 /f a 2 + 295488 a 3 /f vf 2 - 205632 a 3 /f a 2 - 582090 ti /f $ 2 X 

-1157435 ^ 2 $ 2 1 a 2 - 1517760 tt 2 ^ 4 a 2 - 447840 ~f a 2 a 2 

+299000 vr 2 a 2 $ 2 X + 123540 /f $n /5 + 303600 vr /f ac 2 + 51450 ~f a $n 

-126720 vr 2 /f vr 2 - 30643000 tt 3 /?^ - 23040 /f vr 2 a 2 + 78000 ti a 3 $ 2 X 

+311 1840 tt 3 /f W 2 - 361718 p 3 $^ a + 25 p 5 vr $n + 1296000 ti /f /3 2 

-14398260 vr p 2 $ n a a 2 - 57600 p 5 aaW + 781920 ^ a a 2 W 

+89856 a 3 /f a vr + 301945 ]f $ 2 X - 138240 /? 5 (3 a W 

+1886976 J?p a 2 W- 11 19744 pa a 2 + 682560 /f a a (3 

+972820 p 4 a a $n - 3755902 /f a a 2 $u - 266205 ^ 4 $n a ¥ 

+1749108 P A $n a; /3 + 3627158 /f $n a 2 vf - 3766320 /f /3 a 2 $u 

+599880 a 3 $n /f ¥ - 3021720 a 3 $n a /f + 6871680 tt 2 /f ¥ 

-14919600 vr 2 /f $ n /3 - 13552520 tt 2 /f $ n a - 1389792 tt 2 /f a 2 a 

+4357632 tt 2 /f ¥ 2 a + 2915280 tt 2 /f a tt + 13578100 tt 2 /f $n ¥ 

-3718080 n 2 p 4 ap- 60662200 tt 2 /? $ 2 X a - 713952 tt 2 /f vr a a 

-8470120 tt 2 /f $ u fa- 22343840 vr 2 /f $ n a a + 46000 vr 2 a 2 $ n a /? 

-64400 vr 2 a 2 $n /?¥ - 26400 n 2 a $n /3 /f - 105600 vr/f a vf 

+1806460 tt /f $ n a - 38468250 tt /3 $ 2 X a 2 - 1 11360 vr p 4 W 2 a 

+3379920 vr p 4 $n /3 - 474200 7r /? 4 $ n ¥ - 253440 vr /f vr p 

+1268640 np 5 a P - 1647600 Tip 4 a 2 a - 6010420 up 2 <$> 2 n a 

-939744 Tip 3 a 2 a 2 + 1993248 Tip 3 W 2 a 2 + 3021840 Tip 4 a W a 

+51264vr/3 3 ¥aa 2 + 14037554 Tip 3 $ n a ¥ - 14255872 Tip 3 $ n a a 

-15000840 tt p 3 $n a /3 - 592500 ti /f $n a 2 ¥ - 4074624 vr /f a a /3 

+7198848 7r ~p 4 W a P — 14400 ti a 2 $n /? ~p 2 - 16800 vr a 3 $n /?¥ 

+12000 vr a 3 $n a]3 - 9855600 vr 3 /f $n ¥ - 11178800 vr 3 /f $n a 

-929760 vr 3 P 3 ¥ a + 44000 tt 3 a $n ap - 61600 tt 3 a $n /?¥ = . (60) 



3 The case $ n = 

Carminati and McLenaghan [7j used the conditions III and V given in Section |2] to 
prove that Huygens' principle is not satisfied if any of the spin coefficients a, P or ti 
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vanish. We now extend the proof for the case in which a/3ir ^ and $n = 0; i.e., we 
shall prove the following theorem: 

Theorem 6 Let V4 be any space-time which admits a spinor dyad with the properties 

°A-,BB = °aI B b > ( 61 ) 

where I B g is a 2-spinor, and 

^ABCD;EE^^ B ^ C O D O E O E = 0, (62) 

R = 0, ®aba B o A o B = . (63) 

Then the validity of Huygens ' principle for the conformally invariant equation (QJ) 
implies that 

<S>ababO A l B o A T* ± . (64) 

Proof: 

When $n = the quantity N\, given by (|60|). factors in the following form : 



N x := -12/3 (65) 



where 



Px := 12 P (3 + 2ir a + 2 aa + 5 Pa + 6Tir + 2T a, (66) 

p 2 := 1188ppa + 2A0apW + U0]3Wn - 1265^vra 

-2250 J3 p vr + 6830 vr 2 a + 7647 a vr a + 2142 a 2 a 

-690 paa- 10805 vf vr 2 - 11529 a W vr - 3078 W a 2 . (67) 

Let us consider first the case in which p 2 = 0. Applying 5 to (|67p and solving for 
5a, we obtain: 

6a ■= -(690120 p pa 3 - 177100 /f a vr + 2475000 pa vr 2 /3 - 97175 a p 3 a 
-186390 a p 3 p + 1716210 a a 2 P 2 - 1131915 a 2 p 2 W - U70219 pa 3 a 
-3875990 7t 2 P 2 W - 341280 P vr 7j 3 + 5791820 vr 2 a ?J 2 + 9784170 a 3 W vr 
+2573586 a 2 /3 /f + 8639361 a 3 P W + 8903160 it 2 pp 2 - 28683640 vr 3 pa 



-1 1970480 7i a a + 18687060 p z Wa z - 6346350 a d vr a + 9567000 P tt a p 
+2615220 pp n a 2 + 33800 p 3 a W + 61600 p 3 W tt + 90203190 vr 2 a /? W 
- 1 1 19420 a 4 a - 4188240 p 2 If a n + 6302070 a /f vr a 
-24866544 ?J a 2 n a - 46210320 pa tt 2 a + 48333948 a 2 p W tt 
+56154880 7r 3 ^vr + 1705860 a 4 vr - 7513000 tt 3 a a + 11885500 Wn 3 a) j 

((-12a - 22vr + 6c) (1 15?Jc7 + 126/373 - 40^¥ + 783¥a - 1634vra 
+1448vr¥- 921aa)) , (61 
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where the denominator of the expression above, given by 



d A : = (-12a - 22tt + 0) (llhpa + 126/3/3 - 40(3w + 783vra - 1634™ 

+14487T7f - 921aa) , (69) 

is assumed to be non-zero, for now. 

Here Ni, and all equations obtained by comparing different expressions for 5a, 
are polynomials in three complex variables a, (5 and it. One complex variable can be 
eliminated by introducing the following new variables: 

«,:=H, X2 :=L (70) 

TT TT 

In what follows we first prove that the necessary conditions imply that both x\ and 
x 2 are constants. Then, later, we shall prove that this leads to a contradiction. 
In the new variables defined by (J7U|) . the expression (JBTj) assumes the form 

p 2 = -2250x^2 - 1188x2-x 2 a;i - 3078x 2 _ 11529s! - 10805 

+7647xixI+ 6830 + 2142x 2 xlx7- 1265x^x7 + 440 x^ - 690 x^x x x7 
+240xix^ + 2142x^x7 = 0. (71) 

Subtracting ()68j) from (|58jl (with $n = 0), and taking the numerator, gives 

N 2 := -75130000 x7 2 xi -568034312 x^xixT 2 - 162662000 x 2 x7 2 
-263829680 x^x7 2 + 69 8 2 8 00 x7x^x 2 x 1 + 91299060 xf 
-105963000 X0-X2X1 + 328900 x^x7 2 - 3248115 xf - 97977600 x\ 
-12927600 x\ - 278399100 x\ + 37400 xf - 16070400 x\x^x 2 
-277285752 x\ xf x 2 + 6646212 x\ xf x 2 + 78647544 x\xixi 
-48437136 x\xf x 2 + 15461670 x x xlxi 3 - 39517398 x\xf xl 
+ 12081744 xlx\ x? x 2 + 593329572 x?x^ x7 + 4160700 x\ xf x7 
-221403987x 2 x^ 2 x7- 42924720 x7x2" 3 x 2 - 459117172 xT^x 2 
+ 1679968716x7x^x 2 + 112691520 xlx^ 2 x 2 

+ 163297776 x7x2- 2 x 2 xi -413617894 x7x^ 2 xi + 125141836 x7 2 X2~ 2 xi 

+14366310 xlxf - 169509600 x 2 X2"x 2 - 15940800 x\xf x\ 

+114001200 x 2 x^x 2 x7+ 62078400 x 3 x^x 2 x7 

+ 11275200 x\ x^x 2 x7 - 45650088 x 1 xlxf x 2 - 213660 xf x x x 2 

+77697000 x7x 2 X2" 2 x 2 - 12134448 x7x 2 xf x 2 

-165106878 xfx\x 2 - + 66914052 x7 2 x\ xf + 383646000 x\ x7 

+ 136587600x^x7- 1405233644 x 2 ^- 495836310 x 3 ^ 

+78021681 x\ xf - 257647540 x7 x^ 2 - 44763192 x\ xf x x 
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-14850000 x\x x x? - 121050 x^xixT- 11901168 x 2 xi 3 x^ 
-4276800 x\xf x\ - 22290588 xf x\xi + 11929896 xf x\xf 
-3427140 xf x\ x? - 65630268 x\xi+ 13886406 x\ xf 
-939210 x\ xf - 1 1744220 xf xf + 78035680 xf xf + 20400 xf x x 
+ 1533600 x\xf + 146164809 xT 2 xi + 750960 x x xf x 2 xl 
+ 1 79400 ~xf xf X\ - 90396000 x\x^x 2 - 12687610 x x xf xf 
+25076106 xf xi x 2 - 528822552 xf x x x 2 - 3492935 xi 3 x x 
+23644920 xf x 2 - 397830 x^x 2 - 336008520 xf x 2 - 6426000 xf x\ 
-131992500 x\ xf + 18230400 x\xi + 824256 x x xf x\ 
+ 1395480x^x1x2 - 221925x^x1-47574000 x^xT 2 - 166397000 x x 
-42109200 x\ xf - 836819360 xT - 1770562898 x x x 2 ~ + 224037000 x x xl 
+478803300 x\ xl - 351507100 x\ + 999889240 x^xT 

+2115717928 X2"xi xl = . (72) 



We now wish to determine the solutions of the system of algebraic equations 
{p 2 = 0, N 2 = 0}. This may be accomplished in principle using the Grobner basis 
method of Buchberger [T2] as follows. First, we treat the quantities x x , x 2 , xl, x^ as 
independent variables, as view the quantities p 2 , N 2 as polynomials in these indeter- 
minates over the field of rational numbers. (In the subsequent analysis we may use 
the fact that some variables are complex conjugates of each other, but this will not be 
necessary for our immediate purpose.) Then, by computing a Grobner basis for the 
set {p 2 , N 2 } (actually, the ideal < p 2 , N 2 >) with respect to a purely lexicographic 
ordering of terms (see ^2]) we obtain a new set of polynomials with the same solutions 
but in which the variables have been successively eliminated as far as possible. In 
order to speed the computations, we use a special variant of the algorithm jH] which 
combines the nonlinear elimination with factorization of intermediate results. (This 
algorithm is available in the Maple system as the function gsolve.) For the polyno- 
mials {p 2 , N 2 }, the algorithm produces the following components, which collectively 
contain all solutions: 



G 



i •— 



[ -8 + 23 xl, 11 xa" + 8, 66 x x + 125 ] 



(73) 



G 2 := [ 9108 x 2 x^+ 247, -8 + 23 XT, 66xi + 125], 
G 3 := [ 828 x 2 xT + 75 xT + 77, -8 + 23 xl, 66 x x + 125 ] , 
G A := [271xi + 138x 2 xT +517, -8 + 23x1], 
G 5 := [36 x 2 xT + 7 — 5 xl, 6 x x + 11 ] , 



(74) 
(75) 
(76) 
(77) 



G 6 := [671514624 x\ xf + 488374272 x\ xT - 220446720 xTx 2 

+35785728 xf x 2 + 88473600 x 2 xT - 27979776 x x x 2 + 69101568 x 2 
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-167878656 xfxf + 181020672 xT 2 xi - 26599040 xf 

+73852416 xf xl - 132857600 x7 xi + 23168456 x7 - 26978094 x x 

-49722705 - 2204136 x\ + 48043776 x x xi + 111324800 xi - 7645440 xf, 

847872 x7xix 2 - 294912 x 2 xi + 423936 xfx^ - 382720 xf 

-218112 x7xi + 323928x1 - 54450 x x - 169493 + 24576 xi, 

139392 xf x 2 + 202752 xf x 2 + 73728 x 2 xi + 69696 xf xj 

+38456 xf x7 — 54656 xTxi - 33280 xl - 4224 xi + 10432 - 11616 x x xi 

+22544x2"+ 2827 X2" 2 - 11616x^ - 7986 xi 2 Xi, 304128 x x x 2 xi 

+576000 x 2 xi + 66240 xlxi - 59800 xl - 191598 x x - 279575 - 17424 x\ 

+ 198000xi xi + 351960 xi, 192 x x xl - 282 x x - 505 + 280 xl] . (78) 

Using the fact that the pairs (xi, x 2 ) and (x% , xi) are complex conjugates of each 
other, we conclude that the sets G x to G5 provide solutions which are either impossible 
or in which x x and x 2 are constant. In the case of Gq, this is not immediately obvious. 
Its smallest term is: 

192xi xl- 282xi - 505 + 280x7= 0. (79) 

Subtracting (fT9*j) from its complex conjugate we obtain the conclusion that xi is real, 
which implies that it must be constant. It follows that x 2 must be constant as well. 
Let us consider now the case 

Pi = 12x 2 xi + 6 + 2xi + 2x7 + 2xix7 + 5xix 2 = . (80) 

We then use the side relation S x given by (p)2*j) . whose numerator takes the form: 

p 3 ■= -6 x 2 2 xf + 1210 x 2 xi + 1276 x x + 1276 x7 + 360 x 2 x\ xl 

+2901 x 2 xi + 1528 xi x 2 xi + 660 xi xix7 + 660 x 2 x\ + 408 xi x7 2 

+660 x7 2 xi + 408 x7 x\ + 264 x7 2 + 144 x\ xf + 1452 + 264 x\ 

+ 1210 x~i xi + 1444 x7xi + 1528 x 2 x7 3Ji + 660 x 2 x7 xi 

+803x 2 xixix7 + 360xixix7 2 = . (81) 

Applying our nonlinear elimination algorithm as before to p x , p% we obtain the 
following equivalent system of equations: 

s x := 6x 2 xi + 31xix7 + 56xi + x7 + 3 =0, (82) 
s 2 := 72x 2 xix7 + 132xiX 2 + 31xix7 2 + 56x~7xi + x7 2 + 3x7 = , (83) 
S3 := xix 2 — 22xi — 12x7xi = . (84) 

Subtracting (fH2|) from its complex conjugate yields xl = x x . Subtracting (|H4^1 from 
its complex conjugate now gives x 2 = xi = 12x x + 22. Substituting these relations 
back in ()82|) and ()83|) results in a system with no solution. 
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It thus follows that in either of the cases which arise from equation ([65)1 . x\ and 
X2 must necessarily be constant. However, it may be shown (though we postpone the 
details until the following section) that this too leads to a contradiction. 

We must finally consider the case in which the denominator of Dfi, given by ()49j) . is 
zero. Here we shall suppose that $n is not necessarily zero, so that the side relations 
derived in Section El will remain valid in the following section as well. According to 
dUD and (J7UJ), 

d x := 22 + 12xi -xj = 0. (85) 

From (|49|) we obtain 

Ei := -53xi0ii + 242 - 220x 2 + 24x1^ + 30x^011 + 274x2^ 

- 5xi+ 110xiXT + 24 x 2 + 44x1+ 144xi x 2 xi- 108 0n 

— 120 Xix 2 + 68 Xi xi + 176 xi + 60 xTx^xi = , (86) 

where 0n is defined as follows: 

0u := (87) 

7T7T 

Applying 5 to fx, using (jSOJ), (|3~%|) and (f3Hj). and solving for 5a, we get 

5a = 120/3a + 66vra + 220vr/3 + 33a 2 - ~f . (88) 

By applying 5 to d x , now using (pfUj) . (jHJ), (JHJ) and (|4*7|) and solving for D/z, we get 

D//= (47ra + 29^+64aa-370/37r + 327r7r-190/5a + 1147ra + 5^a-68$ii)/20. (89) 

Subtracting given by (jUj), from the complex conjugate of (jHHJ), gives 

E 2 := 20011 - 2x2X^- 22x1 - 24x1x1 + 66x 2 + 35xiX2-22xT 

+35x2" xi +66x2" = 0. (90) 

Applying 5 to (fHUJ) gives 

£ 3 := 1168xi 0ii + 3980x 2 x 2 - 264 + 14520x2 - 3584XTX 2 - 255x2"0n 
-2838 x 2 X2" - 9482 xi xj~ — 24200 xi - 12870 x^xT - 12870 x x xi 
-3784 x 2 - 1210 xf - 4928 xT - 1514 x x x 2 xi + 20 x 2 X2" 2 - 625 X2" 2 xT 
+2046 0ii + 15180 xx x 2 - 7480 x x - 6835 xTx^Xi = . (91) 

Applying nonlinear elimination to dx, fx, Ex, Ex, E2 and £"3 , we find that this system 
has no solution. 

The cases where each of the denominators d 2 , d% and dx, that appeared in the 
preceding equations are zero lead to contradictions, according to [23]. The demon- 
stration of this fact follows the steps described above and will not be presented here 
for brevity. 

Thus, for Huygens' principle to be satisfied on Petrov type III space-times we must 
have $11 7^ 0, and Theorem |3 which states this result in a conformally invariant way, 
is proved. 
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4 The case $ n ^ 



We shall now examine the sole case which remains after the analysis of the previous 
section, namely that in which a/3ir ^ and $n 7^ 0. This, in view of Theorem® will 
complete the proof of Theorem ^ Our approach is related to that of the previous 
section, in that we reduce the problem to an issue of solvability of a purely algebraic 
system of equations. 

We first observe that, in addition to the algebraic equations given by (JH2*|) and 
(jnOJ), an extra independent equation may be obtained by applying the NP operator 5 
to (|52p. All of the Pfaflians which result are known explicitly, and may be replaced 
using the expressions found in Section to obtain a (very large) expression in the 
complex variables a, /3, and 7r and the real quantity $n. Upon transforming variables 
according to (JTOj) and (|S7jl . we obtain a complex quantity in the new variables x%, 
X2, xi, x~2, and <pi\. (This polynomial contains 408 terms of maximum total degree 
9.) Together with the equations which similarly follow from (J52j) and (|6U|) . we have 
in effect a system of five equations in five real variables. Let us denote the set of 
polynomials which arise in this system (i.e. when the equations are written with a 
right hand side of 0) by F . 

It must be mentioned that the approach of the previous section, namely computing 
the solutions by explicit elimination, is impossible in the present case due to the 
intrinsic computational complexity of nonlinear elimination and the high degree of 
our polynomials. It is possible and will suffice, however, to bound the number of 
solutions using the following result due to Buchberger [T2] : 

Theorem 7 Let G be a Grobner basis for < F > (the polynomial ideal generated 
by F) with respect to a given ordering of terms, and let H denote the set of leading 
terms of the elements of G with respect to the chosen term ordering. Then the system 
of equations corresponding to F has finitely many solutions if and only if for every 
indeterminate x in F there is a natural number m such that x m G H. 

The key to using this result is that we may use an ordering of terms based on total 
degree (i.e. a non-elimination ordering) for which the computational complexity of 
Buchberger's algorithm for Grobner bases is much lower. Unfortunately, even in this 
setting a Grobner basis for F cannot easily be computed due to the extreme size of 
intermediate results produced by the algorithm. 

It would be highly desirable to apply modular homomorphisms in the manner 
used in algorithms for factorization (e.g. so-called Chinese remainder, or Hensel 
algorithms J2]) i n the present situation. This is not currently possible due to a 
number of unresolved problems with the approach. Nonetheless, it provides a useful 
probabilistic experimental approach: treat the elements of F as polynomials over a 
prime field Z p (rather than the rationals), where p is of modest size, and compute 
the Grobner basis of F modulo p over Z p . For a single prime, it is possible that the 
result so obtained may have no useful relationship with the Grobner basis of F over 
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the rationals. However, if the basis polynomials computed using a large number of 
different primes all exhibit identical monomial structures, it is extremely likely that 
they each represent a distinct homomorphic image of the true Grobner basis of F. 
The question of accurately computing the probability of success for a specific series of 
primes remains an open problem. However, the individual prime field computations 
are comparatively easy since (unlike the rational case) no single coefficient may be 
larger than the chosen prime. This provides an experimental "sampling" method 
which gives clues on how best to compute the true result, and what that result will 
likely be. 

We must also consider that if we were able to compute a Grobner basis for F 
over the rationals, we would derive information on all solutions of the corresponding 
system including those which were examined in the previous section (i.e. for which 
$n = 0). It is possible to exclude those solutions entirely by adding an additional 
constraint and variable, 

(p n z -1 = (92) 

to our equations to produce the augmented system F. Still, only an actual computa- 
tion reveals whether this improves or worsens the tractability of the problem. In our 
case, a large number (a few thousand) prime field "sample" computations (done using 
the GB package of Faugere 11J, which is far more efficient than the general-purpose 
Maple system) all suggested that the addition of equation (J9*2~j) made the Grobner 
basis calculation much more efficient. More importantly, once the solutions examined 
in the previous section were in effect discarded, only a finite number remained when 
Theorem [7| is taken into account. With this in mind, it was possible (and worthwhile) 
to compute the true Grobner basis of F over the rationals in the indeterminates x\, 
X2, xT, ~x~2, Tn, z using a total degree ordering of terms. Since this basis contains 
polynomials with leading terms 

Xi6, #25, xT5, x~2~5, 0ii, z5, (93) 

we may conclude that there are only finitely many solutions for which 0n, and hence 
$n as well, is nonzero. (For this last computation the latest and most efficient version 
of Faugere's GB package, known as FGB, was required.) It follows that x\, x 2 , xT, 
T2, 0n must be constants; it remains only to show that this yields a contradiction. 

Since 0n must be constant (including the case in which $ n = 0) it follows from 
(|7Uj). (|57jl that the quantities ttW, (3(3 and a/ (3 are all constant as well. From the 
equation S(/3(3) = we obtain, in the variables x±, x 2 , 0n given by (f70|l . ljB7j) . the 
side relation 

(7x7. + 25x 2 + 2)0n + 374x2X2 + 199x7x2X2 — 5x2X2 
+60xixTx 2 + 60x~T 2 x~2 + 110xix 2 + 110xiX 2 + 68x1x7 

+24xixT 2 + 24xl 2 + 44xi + 116x1 + 132 = . (94) 
Next, from <5(x2) = we obtain the Pfaffian 

Sw = -W(a + P) . (95) 
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Using this, along with the previously determined Pfaffians, we then obtain from 
5(nW) = another side relation; on subtracting this result from (}9*4*|) (and ignoring 
the possibility that d% = 0, which has already been considered) we obtain 

x 2 (xi + x^ + 4) + 0H = 0. (96) 

Finally, from 5 (a/ f3) = we obtain 

+ 5x^ + 2) + 9x5 = 0. (97) 

The collection of polynomials given by (JMJ), (}96|) . (|9"?|) . (|52|) and their complex con- 
jugates has a Grobner basis (computed easily using Maple) containing only the poly- 
nomial 1. This is equivalent to showing that there exists a combination of these 
polynomials which equals 1, and hence that they cannot vanish simultaneously (see 
[12]); i.e., the associated system of equations has no solutions. This completes the 
proof. 

5 Conclusion 

In completing the proof of Theorem Q we have fully solved Hadamard's problem for 
the scalar wave equation in the case of Petrov type III space-times. Essential to our 
proof were use of the six-index necessary condition obtained by Rinke and Wiinsch 
|221) an d separate analyses (and different ideal-theoretic tools) for the cases $n = 
and $n 7^ 0. To complete the proof of the conjecture stated in the Introduction it 
remains to consider the space-times of Petrov types I and II. A partial result for type 
II has been obtained by Carminati, Czapor, McLenaghan and Williams However, 
it is not yet clear whether the complicated equations which arise from conditions III, 
V, and VII can be solved by the method used in the present paper. 
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